Genomes of multicellular algal sisters to land plants illuminate signaling network evolution.

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Aethionema arabicum dimorphic seed trait resetting during transition to seedlings.

The transition from germinating seeds to emerging seedlings is one of the most vulnerable plant life cycle stages. Heteromorphic diaspores (seed and fruit dispersal units) are an adaptive bet-hedging strategy to cope with spatiotemporally variable environments. While the roles and mechanisms of seedling traits have been studied in monomorphic species, which produce one type of diaspore, very little is known about seedlings in heteromorphic species. Using the dimorphic diaspore model Aethionema arabicum (Brassicaceae), we identified contrasting mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained from IND fruits by pericarp (fruit coat) removal. What follows the completion of germination is the pre-emergence seedling growth phase, which we investigated by comparative growth assays of early seedlings derived from the M+ seeds, bare M- seeds, and IND fruits. The dimorphic seedlings derived from M+ and M- seeds did not differ in their responses to ambient temperature and water potential. The phenotype of seedlings derived from IND fruits differed in that they had bent hypocotyls and their shoot and root growth was slower, but the biomechanical hypocotyl properties of 15-day-old seedlings did not differ between seedlings derived from germinated M+ seeds, M- seeds, or IND fruits. Comparison of the transcriptomes of the natural dimorphic diaspores, M+ seeds and IND fruits, identified 2,682 differentially expressed genes (DEGs) during late germination. During the subsequent 3 days of seedling pre-emergence growth, the number of DEGs was reduced 10-fold to 277 root DEGs and 16-fold to 164 shoot DEGs. Among the DEGs in early seedlings were hormonal regulators, in particular for auxin, ethylene, and gibberellins. Furthermore, DEGs were identified for water and ion transporters, nitrate transporter and assimilation enzymes, and cell wall remodeling protein genes encoding enzymes targeting xyloglucan and pectin. We conclude that the transcriptomes of seedlings derived from the dimorphic diaspores, M+ seeds and IND fruits, undergo transcriptional resetting during the post-germination pre-emergence growth transition phase from germinated diaspores to growing seedlings.

The dimorphic diaspore model Aethionema arabicum (Brassicaceae): Distinct molecular and morphological control of responses to parental and germination temperatures.

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting ABA sensitivity. This involved expression of morph-specific transcription factors, hypoxia response and cell wall-remodeling genes, as well as altered abscisic acid (ABA) metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

Near telomere-to-telomere genome of the model plant Physcomitrium patens.

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Economic benefits of herd genotyping and using sexed semen for pure and beef-on-dairy breeding in dairy herds.

The cost benefits of herd genotyping and the benefits of using sexed semen have been affected by recent improvements in sexing technologies, incorporation of direct health traits in the German total merit index for Holstein cattle, deteriorating prices for purebred heifer calves and bull calves, and introduction of herd genotyping programs. Inseminating genetically superior dams with female-sexed Holstein semen increases the mean breeding value of heifer calves and can produce more Holstein heifer calves than are needed for replacement. This provides an opportunity to increase the selection response in health and production traits at the farm level. A deterministic model is introduced that predicts the increase or decrease in net profit when a farmer takes part in a herd genotyping program and follows a certain insemination strategy. The types of semen allocated to cows and heifers may be sexed or unsexed and Holstein or beef breed. Genetically superior heifers and cows are inseminated with female-sexed Holstein semen, intermediate dams with unsexed Holstein semen, and genetically inferior dams with unsexed or male-sexed beef breed semen. In general, participating in a herd genotyping program is beneficial for German Holstein breeders. The optimum proportions of cows and heifers that should be inseminated with a certain type of semen are sensitive to farm-specific peculiarities. A small price difference between crossbred bull calves and crossbred heifer calves often makes the use of male-sexed beef breed semen uneconomic. Under the conditions considered, it was found to be advantageous to inseminate approximately 50% of heifers and 10% of cows with the highest genetic merit with female-sexed Holstein semen. The optimum proportion of cows that should be inseminated with unsexed beef breed semen was found to be approximately 40%. In a herd with a low replacement rate, the selected heifers can exhibit their genetic superiority over a longer period of time, and a larger proportion of cows can be inseminated with beef breed semen. Participation in a herd genotyping program is, therefore, particularly beneficial for herds with low replacement rates.

Genomic imprints of unparalleled growth.

Chlorella ohadii was isolated from desert biological soil crusts, one of the harshest habitats on Earth, and is emerging as an exciting new green model for studying growth, photosynthesis and metabolism under a wide range of conditions. Here, we compared the genome of C. ohadii, the fastest growing alga on record, to that of other green algae, to reveal the genomic imprints empowering its unparalleled growth rate and resistance to various stressors, including extreme illumination. This included the genome of its close relative, but slower growing and photodamage sensitive, C. sorokiniana UTEX 1663. A larger number of ribosome-encoding genes, high intron abundance, increased codon bias and unique genes potentially involved in metabolic flexibility and resistance to photodamage are all consistent with the faster growth of C. ohadii. Some of these characteristics highlight general trends in Chlorophyta and Chlorella spp. evolution, and others open new broad avenues for mechanistic exploration of their relationship with growth. This work entails a unique case study for the genomic adaptations and costs of exceptionally fast growth and sheds light on the genomic signatures of fast growth in photosynthetic cells. It also provides an important resource for future studies leveraging the unique properties of C. ohadii for photosynthesis and stress response research alongside their utilization for synthetic biology and biotechnology aims.

Proteome plasticity during Physcomitrium patens spore germination - from the desiccated phase to heterotrophic growth and reconstitution of photoautotrophy.

The establishment of moss spores is considered a milestone in plant evolution. They harbor protein networks underpinning desiccation tolerance and accumulation of storage compounds that can be found already in algae and that are also utilized in seeds and pollen. Furthermore, germinating spores must produce proteins that drive the transition through heterotrophic growth to the autotrophic plant. To get insight into the plasticity of this proteome, we investigated it at five timepoints of moss (Physcomitrium patens) spore germination and in protonemata and gametophores. The comparison to previously published Arabidopsis proteome data of seedling establishment showed that not only the proteomes of spores and seeds are functionally related, but also the proteomes of germinating spores and young seedlings. We observed similarities with regard to desiccation tolerance, lipid droplet proteome composition, control of dormancy, and ß-oxidation and the glyoxylate cycle. However, there were also striking differences. For example, spores lacked any obvious storage proteins. Furthermore, we did not detect homologs to the main triacylglycerol lipase in Arabidopsis seeds, SUGAR DEPENDENT1. Instead, we discovered a triacylglycerol lipase of the oil body lipase family and a lipoxygenase as being the overall most abundant proteins in spores. This finding indicates an alternative pathway for triacylglycerol degradation via oxylipin intermediates in the moss. The comparison of spores to Nicotiana tabacum pollen indicated similarities for example in regards to resistance to desiccation and hypoxia, but the overall developmental pattern did not align as in the case of seedling establishment and spore germination.

A DELAY OF GERMINATION 1 (DOG1)-like protein regulates spore germination in the moss Physcomitrium patens.

DELAY OF GERMINATION 1 is a key regulator of dormancy in flowering plants before seed germination. Bryophytes develop haploid spores with an analogous function to seeds. Here, we investigate whether DOG1 function during germination is conserved between bryophytes and flowering plants and analyse the underlying mechanism of DOG1 action in the moss Physcomitrium patens. Phylogenetic and in silico expression analyses were performed to identify and characterise DOG1 domain-containing genes in P. patens. Germination assays were performed to characterise a Ppdog1-like1 mutant, and replacement with AtDOG1 was carried out. Yeast two-hybrid assays were used to test the interaction of the PpDOG1-like protein with DELLA proteins from P. patens and A. thaliana. P. patens possesses nine DOG1 domain-containing genes. The DOG1-like protein PpDOG1-L1 (Pp3c3_9650) interacts with PpDELLAa and PpDELLAb and the A. thaliana DELLA protein AtRGA in yeast. Protein truncations revealed the DOG1 domain as necessary and sufficient for interaction with PpDELLA proteins. Spores of Ppdog1-l1 mutant germinate faster than wild type, but replacement with AtDOG1 reverses this effect. Our data demonstrate a role for the PpDOG1-LIKE1 protein in moss spore germination, possibly alongside PpDELLAs. This suggests a conserved DOG1 domain function in germination, albeit with differential adaptation of regulatory networks in seed and spore germination.

How economic weights translate into genetic and phenotypic progress, and vice versa.

BACKGROUND: This paper highlights the relationships between economic weights, genetic progress, and phenotypic progress in genomic breeding programs that aim at generating genetic progress in complex, i.e., multi-trait, breeding objectives via a combination of estimated breeding values for different trait complexes. RESULTS: Based on classical selection index theory in combination with quantitative genetic models, we provide a methodological framework for calculating expected genetic and phenotypic progress for all components of a complex breeding objective. We further provide an approach to study the sensitivity of the system to modifications, e.g. to changes in the economic weights. We propose a novel approach to derive the covariance structure of the stochastic errors of estimated breeding values from the observed correlations of estimated breeding values. We define 'realized economic weights' as those weights that would coincide with the observed composition of the genetic trend and show, how they can be calculated. The suggested methodology is illustrated with an index that aims at achieving a breeding goal composed of six trait complexes, that was applied in German Holstein cattle breeding until 2021. CONCLUSIONS: Based on the presented results, the main conclusions are (i) the composition of the observed genetic progress matches the expectations well, with predictions being slightly better when the covariance of estimation errors is taken into account; (ii) the composition of the expected phenotypic trend deviates significantly from the expected genetic trend due to the differences in trait heritabilities; and (iii) the realized economic weights derived from the observed genetic trend deviate substantially from the predefined ones, in one case even with a reversed sign. Further results highlight the implications of the change to a modified breeding goal based on the example of a new index comprising eight, partly new, trait complexes, which is used since 2021 in the German Holstein breeding program. The proposed framework and the analytical tools and software provided will be useful to define more rational and generally accepted breeding objectives in the future.

Chromosome-level genomes of multicellular algal sisters to land plants illuminate signaling network evolution.

The filamentous and unicellular algae of the class Zygnematophyceae are the closest algal relatives of land plants. Inferring the properties of the last common ancestor shared by these algae and land plants allows us to identify decisive traits that enabled the conquest of land by plants. We sequenced four genomes of filamentous Zygnematophyceae (three strains of Zygnema circumcarinatum and one strain of Z. cylindricum) and generated chromosome-scale assemblies for all strains of the emerging model system Z. circumcarinatum. Comparative genomic analyses reveal expanded genes for signaling cascades, environmental response, and intracellular trafficking that we associate with multicellularity. Gene family analyses suggest that Zygnematophyceae share all the major enzymes with land plants for cell wall polysaccharide synthesis, degradation, and modifications; most of the enzymes for cell wall innovations, especially for polysaccharide backbone synthesis, were gained more than 700 million years ago. In Zygnematophyceae, these enzyme families expanded, forming co-expressed modules. Transcriptomic profiling of over 19 growth conditions combined with co-expression network analyses uncover cohorts of genes that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

DELLA proteins regulate spore germination and reproductive development in Physcomitrium patens.

Proteins of the DELLA family integrate environmental signals to regulate growth and development throughout the plant kingdom. Plants expressing non-degradable DELLA proteins underpinned the development of high-yielding 'Green Revolution' dwarf crop varieties in the 1960s. In vascular plants, DELLAs are regulated by gibberellins, diterpenoid plant hormones. How DELLA protein function has changed during land plant evolution is not fully understood. We have examined the function and interactions of DELLA proteins in the moss Physcomitrium (Physcomitrella) patens, in the sister group of vascular plants (Bryophytes). PpDELLAs do not undergo the same regulation as flowering plant DELLAs. PpDELLAs are not degraded by diterpenes, do not interact with GID1 gibberellin receptor proteins and do not participate in responses to abiotic stress. PpDELLAs do share a function with vascular plant DELLAs during reproductive development. PpDELLAs also regulate spore germination. PpDELLAs interact with moss-specific photoreceptors although a function for PpDELLAs in light responses was not detected. PpDELLAs likely act as 'hubs' for transcriptional regulation similarly to their homologues across the plant kingdom. Taken together, these data demonstrate that PpDELLA proteins share some biological functions with DELLAs in flowering plants, but other DELLA functions and regulation evolved independently in both plant lineages.

Co-action of COP1, SPA and cryptochrome in light signal transduction and photomorphogenesis of the moss Physcomitrium patens.

The Arabidopsis COP1/SPA ubiquitin ligase suppresses photomorphogenesis in darkness. In the light, photoreceptors inactivate COP1/SPA to allow a light response. While SPA genes are specific to the green lineage, COP1 also exists in humans. This raises the question of when in evolution plant COP1 acquired the need for SPA accessory proteins. We addressed this question by generating Physcomitrium Ppcop1 mutants and comparing their visible and molecular phenotypes with those of Physcomitrium Ppspa mutants. The phenotype of Ppcop1 nonuple mutants resembles that of Ppspa mutants. Most importantly, both mutants produce green chloroplasts in complete darkness. They also exhibit dwarfed gametophores, disturbed branching of protonemata and absent gravitropism. RNA-sequencing analysis indicates that both mutants undergo weak constitutive light signaling in darkness. PpCOP1 and PpSPA proteins form a complex and they interact via their WD repeat domains with the VP motif of the cryptochrome CCE domain in a blue light-dependent manner. This resembles the interaction of Arabidopsis SPA proteins with Arabidopsis CRY1, and is different from that with Arabidopsis CRY2. Taken together, the data indicate that PpCOP1 and PpSPA act together to regulate growth and development of Physcomitrium. However, in contrast to their Arabidopsis orthologs, PpCOP1 and PpSPA proteins execute only partial suppression of light signaling in darkness. Hence, additional repressors may exist that contribute to the repression of a light response in dark-exposed Physcomitrium.

Genetic survey in smallholder dairy production system of Southern province of Zambia.

Milk production among the smallholder dairy farmers in Zambia is reported to be low despite improvements in milk collection infrastructure and unmet demand by the populace. This study was conducted to characterize cattle breeds and the breeding strategies of the smallholder dairy production system in the Southern province of Zambia. The study was conducted using questionnaires to obtain responses from identified respondents in six districts considered the main dairy cattle producing areas. One hundred and twenty questionnaires were administered and available for descriptive statistical analysis. The cattle genetic groups included local (Tonga and Angoni) breeds and their crosses (40.0%); beef (Boran and Brahman) breeds and their crosses (23.3%); Friesian breed and its crosses (28.8%); Jersey breed and its crosses (4.9%); and Fleckvieh breed (3.0%). Calving rate was found to be to be between 36% and 50% for the dairy herds. Compared to the other genetic groups, it was noted that Friesian and its crosses had higher milk production (3.55 L vs. 1.4 L). The low milk yield reported could be attributed to the period of study, which coincided with the long dry season between April and October. The Friesian and its crosses recorded earlier age at first calving (2.65 vs. 2.8 years), shorter longevity (10.5 vs. 12.25 years), and shorter lactation length (7.5 vs. 9.25 months). Jersey and its crosses, on the other hand, had an above average performance for all indices of economic importance such as milk production, lactation length, age at first calving, longevity, and number of calves produced. The results indicated natural mating was practised by 85% of the farmers. Farmer-preferred traits include size and conformation (23%); performance and colour (16% and 15%, respectively) in the choice of a bull. Selection by farmers were thus based on simple observation without pedigree or performance-based genetic evaluation. It was apparent that the smallholder dairy production system is in a shifting trend to transform the mainly local and beef breed animals into dairy herds through crossing with exotic dairy breeds. The production system is, however, faced with the challenges of water and feed supply to meet nutritional requirements, and high disease burden. Jersey breed was found to be a logical choice for the resource poor smallholder dairy farmers. The development and management of the smallholder dairy breeding schemes should be all-inclusive and directed at the prevalent production systems with the aim of also improving the feeding and management practices.

Establishment and optimization of a new model organism to study early land plant evolution: Germination, cultivation and oospore variation of Chara braunii Gmelin, 1826.

For studying land plant evolution, the establishment and optimization of model organisms representing streptophytic algae, sister to land plants, is essential. Long-term cultivation experiments with Chara braunii S276 were performed over 8 years, since 4 years (Nov. 2018) under constant conditions. Additionally, short-term experiments for optimization of culture conditions were performed with three strains of C. braunii (S276, NIES-1604 and Lausiger Teiche, LaT-2708). Germination success after application of sterilization agents, addition of gibberellic acid and under different incubation conditions with respect to pre-treatment, irradiance regime and substrate was investigated in order to develop protocols for generative cultivation of at least unialgal cultures. The resulting cultivation protocols for C. braunii S276, allowing maintenance of vegetative as well as generative cultures are presented in detail, including protocols for germination induction and growth of sterilized and unsterilized oospores.

Comparative transcriptomics identifies candidate genes involved in the evolutionary transition from dehiscent to indehiscent fruits in Lepidium (Brassicaceae).

BACKGROUND: Fruits are the seed-bearing structures of flowering plants and are highly diverse in terms of morphology, texture and maturation. Dehiscent fruits split open upon maturation to discharge their seeds while indehiscent fruits are dispersed as a whole. Indehiscent fruits evolved from dehiscent fruits several times independently in the crucifer family (Brassicaceae). The fruits of Lepidium appelianum, for example, are indehiscent while the fruits of the closely related L. campestre are dehiscent. Here, we investigate the molecular and genetic mechanisms underlying the evolutionary transition from dehiscent to indehiscent fruits using these two Lepidium species as model system. RESULTS: We have sequenced the transcriptomes and small RNAs of floral buds, flowers and fruits of L. appelianum and L. campestre and analyzed differentially expressed genes (DEGs) and differently differentially expressed genes (DDEGs). DEGs are genes that show significantly different transcript levels in the same structures (buds, flowers and fruits) in different species, or in different structures in the same species. DDEGs are genes for which the change in expression level between two structures is significantly different in one species than in the other. Comparing the two species, the highest number of DEGs was found in flowers, followed by fruits and floral buds while the highest number of DDEGs was found in fruits versus flowers followed by flowers versus floral buds. Several gene ontology terms related to cell wall synthesis and degradation were overrepresented in different sets of DEGs highlighting the importance of these processes for fruit opening. Furthermore, the fruit valve identity genes FRUITFULL and YABBY3 were among the DEGs identified. Finally, the microRNA miR166 as well as the TCP transcription factors BRANCHED1 (BRC1) and TCP FAMILY TRANSCRIPTION FACTOR 4 (TCP4) were found to be DDEGs. CONCLUSIONS: Our study reveals differences in gene expression between dehiscent and indehiscent fruits and uncovers miR166, BRC1 and TCP4 as candidate genes for the evolutionary transition from dehiscent to indehiscent fruits in Lepidium.

The greening ashore.

More than half a billion years ago a streptophyte algal lineage began terraforming the terrestrial habitat and the Earth's atmosphere. This pioneering step enabled the subsequent evolution of all complex life on land, and the past decade has uncovered that many traits, both morphological and genetic, once thought to be unique to land plants, are conserved across some streptophyte algae. They provided the common ancestor of land plants with a repertoire of genes, of which many were adapted to overcome the new biotic and abiotic challenges. Exploring these molecular adaptations in non-tracheophyte species may help us to better prepare all green life, including our crops, for the challenges precipitated by the climate change of the Anthropocene because the challenges mostly differ by the speed with which they are now being met.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes.

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Delayed differentiation of epidermal cells walls can underlie pedomorphosis in plants: the case of pedomorphic petals in the hummingbird-pollinated Caiophora hibiscifolia (Loasaceae, subfam. Loasoideae) species.

BACKGROUND: Understanding the relationship between macroevolutionary diversity and variation in organism development is an important goal of evolutionary biology. Variation in the morphology of several plant and animal lineages is attributed to pedomorphosis, a case of heterochrony, where an ancestral juvenile shape is retained in an adult descendant. Pedomorphosis facilitated morphological adaptation in different plant lineages, but its cellular and molecular basis needs further exploration. Plant development differs from animal development in that cells are enclosed by cell walls and do not migrate. Moreover, in many plant lineages, the differentiated epidermis of leaves, and leaf-derived structures, such as petals, limits organ growth. We, therefore, proposed that pedomorphosis in leaves, and in leaf-derived structures, results from delayed differentiation of epidermal cells with respect to reproductive maturity. This idea was explored for petal evolution, given the importance of corolla morphology for angiosperm reproductive success. RESULTS: By comparing cell morphology and transcriptional profiles between 5 mm flower buds and mature flowers of an entomophile and an ornitophile Loasoideae species (a lineage that experienced transitions from bee- to hummingbird-pollination), we show that evolution of pedomorphic petals of the ornithophile species likely involved delayed differentiation of epidermal cells with respect to flower maturity. We also found that developmental mechanisms other than pedomorphosis might have contributed to evolution of corolla morphology. CONCLUSIONS: Our results highlight a need for considering alternatives to the flower-centric perspective when studying the origin of variation in flower morphology, as this can be generated by developmental processes that are also shared with leaves.

The EpiDiverse Plant Epigenome-Wide Association Studies (EWAS) Pipeline.

Bisulfite sequencing is a widely used technique for determining DNA methylation and its relationship with epigenetics, genetics, and environmental parameters. Various techniques were implemented for epigenome-wide association studies (EWAS) to reveal meaningful associations; however, there are only very few plant studies available to date. Here, we developed the EpiDiverse EWAS pipeline and tested it using two plant datasets, from P. abies (Norway spruce) and Q. lobata (valley oak). Hence, we present an EWAS implementation tested for non-model plant species and describe its use.